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Abstract 

The concept of the potential acting between charged particles is reviewed within the 
framework of quantum electrodynamics. The one-photon and two-photon exchange poten¬ 
tials Vi 7 and V 27 are discussed, with attention to the following features: the dependence 
of V 27 on the choice of Ti 7 ; the distinction between FGI (Feynman gauge inspired) and 
CGI (Coulomb gauge inspired) potentials, and the dependence of the large-r behavior on 
these choices. A summary of available results for two spin-0 particles, for a spin-0 and a 
spin 1/2 particle and for two spin-1/2 particles is given. The application of such ’’scattering 
potentials” to the calculation of energy levels of bound states is discussed. 
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1 Introduction and historical review 

It is a great pleasure for me to be able to contribute to this volume in honor of my good friend, 
Ingvar Lindgren. A large part of Ingvar’s illustrious career has been devoted to work in theoretical 
atomic physics. Since the structure and properties of atoms are largely determined by the nature of 
the interaction between its electrons, I thought this would be a nice occasion to present a discussion 
of this topic. I will focus on aspects which are of especial importance for the relativistic theory 
of atomic structure, to which Ingvar has made so many contributions. Because he has had an 
amazingly large number of graduate students (over 60!), there is a good chance that some of them 
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will find something here which Ingvar would think they ought to know. (Note: The exclamation 
point after 60 indicates astonishment, not “factorial” - even Ingvar has his limitations.) 

The question posed in the title may be short, but it does not have a snappy answer. If one 
asks “What is the force between planets?”, everybody knows that it is given with great accuracy 
by Newton’s law of universal gravitation. But for electrons, the story is much more complicated 
and indeed the question itself must be recast. There are at least three reasons for this: (i) physical 
systems involving electrons must be described by quantum mechanics (QM), (ii) in many such 
systems, the relative velocities are either not very small compared to c and/or high accuracy is of 
interest, so special relativity must be used, and (iii) electrons have spin. 

With regard to (i), of course in nonrelativistic classical mechanics the concept of force and, in 
particular, the concept of the force between particles is well-dehned; if the force is conservative, 
so is the concept of two-body potential U(l,2). Thus, in this context the answer is the familiar 
one. The force is the Coulomb force, known for 200 years, and for two point particles, “1” and 
“2”, the associated potential is just the Coulomb potential 17c'(l,2), 

f/c(l, 2 ) =/cia/r, [/ci 2 = eiCa/dvr] ( 1 ) 

with r = |ri — r 2 |; for two electrons, ki 2 a = e^/47r ~ 1/137. However, in quantum theory the 
concept of force plays a secondary role and the question should really be: “What is the potential 
acting between two electrons?” In the context of nonrelativistic quantum mechanics (NRQM) 
and a conhguration space description of the system in question, the answer is again Uc{l, 2), now 
reinterpreted as a multiplicative operator acting on a many-body Schrodinger wave function. 

With regard to (ii), when the effects of the hnite speed of propagation of light are taken into 
account, things get more complicated, both on conceptual and technical grounds. Now even in 
classical electrodynamics the question becomes murky: While the meaning of the force on “1” is 
sharp, via the Lorentz force law, the force exerted by “2” on “1” depends on the previous history 
of “2” (and vice versa ). The question arises: To what extent can one describe the force between 
“1” and “2” in terms of properties of their motion at a given instant? This question appears to 
have been hrst addressed by C.E. Darwin [1], the grandson of Charles Darwin. Darwin showed 
that, to order n^/c^, the effects of retardation can be taken into account by adding to the free 
Lagrangian for two point particles not only the Coulomb term -Uc but also a term -Uu: with 

Ud = -l/ 2 (vi ■ V2 -h vi ■ fv2 ■ f)Uc/c^ ( 2 ) 

from which the corresponding force may be calculated. On passing to the Hamiltonian formalism 
and introducing canonical momenta one gets a momentum-dependent interaction term of the form 

Ud = -l/2(pi ■ P2 + Pi ■ fp2 ■ f)Uc/mim2C^. (3) 

This must have a counterpart in quantum theory. But how is one to order the factors, if Pi is 
replaced by p°^? In whatever way this question is answered, it becomes clear that at this level 
any effective potential between electrons will not be strictly local. 

Finally, with regard to (iii), after the discovery of the spin-1/2 character of the electron and 
the associated magnetic moment, it was obvious that spin-dependent terms must be added to Ud, 
corresponding to the interaction of the magnetic held produced by “1” with the magnetic moment 
of “2” and vice versa . This led to the now familiar spin-other-orbit and spin-spin potentials and 
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gave rise to an operator which could describe the effect of electron-electron interaction on 
atomic fine structure correctly to order a^Ry (apart from a missing contact term). So one could 
say with some justification that, to order and the operator Uq + describes the force 

between electrons. 

After the development of the Dirac equation one might have guessed that, within a framework 
in which the state of a many-electron system is described by a multi-Dirac spinor, the velocity 
factors Vi in (2) should simply be replaced by their formal counterparts in Dirac theory, viz . cai. 
This yields the operator 


Ub{1, 2) = -l/2(ai ■a 2 + ai ■ TO 2 ■ f)Uc, (4) 

known as the Breit operator. It was first obtained by Breit [2], but not in the way I have described. 
After the development of QED by Dirac, Breit studied the level shift in the helium atom arising 
from the exchange of a transverse photon between the two electrons, within a formally relativistic 
(but actually unsound) framework. He concluded that, to a good approximation, this shift is 
given by the expectation value of Ub with a Dirac-type wave function which includes the effect of 
Uc- From this one would expect that an effective interaction between electrons which takes into 
account items (i), (ii), and (iii) is given by the Coulomb-Breit potential. 


Dcb(1,2) = Dc(1,2) + Db(1,2). (5) 

However, although accepted for many years, this expectation turns out to be misleading at best 
and wrong at worst. 

In the following sections I will consider the question of the electron-electron interaction from 
the perspective of modern QED. It will be useful to consider both spin-1/2 and spin-0 particles. 
Sec. H deals with the relatively simple case of one-photon exchange potentials, in the context of 
the scattering of two point particles. However, 200 years after Coulomb one ought to have some 
idea about the extent to which the exchange of two photons can also be represented by an effective 
potential. This is discussed Sec. HI. In Sec. IV I turn to the question of the connection between 
these potentials and the properties of bound states. 

2 Scattering potentials: One-photon-exchange 

Within the framework of perturbative quantum field theory, there is a sharp contrast between the 
beautiful methods available for the calculation of collision amplitudes and those used in practice 
for the calculation of the properties of composite systems or “bound states,” especially in the case 
of a gauge theory such as QED. As an example, consider the scattering of two particles, “1” and 

1 + 2 ^ 1 + 2 , ( 6 ) 

with initial and final four-momenta pi, p 2 and p[, p' 2 , respectively. In a Lorentz-invariant theory the 
transition amplitude has the form tfi = NfMNi, where Nf = (dEjE^)”^/^ and W = {4EiE2)~^^‘^ 
are the usual kinematic factors, and 

M = M{s,t) [s = P^,t = Q^] (7) 
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is the invariant Feynman amplitude, with s and t the invariant squares of the total four-momentum 
P = Pi + P 2 and four-momentum transfer Q = Pi — p'l- Let us hrst consider the effects arising 
from one-photon exchange between two distinct point-like spin- 1/2 particles, such as an e~ and a 

2.1 Two spin-1/2 particles 

Using Feynman gauge for the photon propagator one hnds that the lowest-order contribution to 
M, regarded as a power series in 6162 , is given by 

= u\u^F^‘^^UiU2 ( 8 ) 

with the Ui denoting lepton spinors normalized to 2 mj and 

= 616271 ■ 72/t. ( 9 ) 

In the c.m. system, Q (0, q), with q the three-momentum transfer, q = p — On taking the 
Fourier transform of F/® with a factor exp{—ic{ ■ r), one sees that the corresponding transition 
amplitude can be obtained by taking the matrix element between plane wave spinors of the 
operator 

Ucg{^i 2) = ki 2 {l — ai ■ 0 . 2 )!^ = Uc(l, 2) -|- Ug(1, 2) (10) 

with Uq = —/i;i 2 (q;i ■ Q; 2 /r), the so-called Gaunt potential [3]. 

It is therefore tempting to say that, to leading order in 6162 , but regardless of the relative speed, 
the effective interaction between, “1” and “2” is given by (10). However, suppose that instead 
one uses the Coulomb gauge for the photon propagator. Then one finds that (9) is replaced by 

Ft ^ = ei62[7i72/Q^ + (71 ■ 72 — 7i ■ QI 2 ■ QIQ^)IQ^\- (11) 

( 2 ) 

Fourier transformation of F^.^ yields 


Ucb{^-,‘ 2) = Uc + Ub, ( 12 ) 

where Us is the operator 

Ub =-kuiai ■ a2 + air ■ a2 ■ r)/ 2 T (13) 

originally obtained by Breit in the context of a bound state problem, as mentioned above. By 
construction, the Coulomb-Gaunt potential Ucg and the Coulomb-Breit potential Ucb have the 
same matrix element between product plane wave spinors, provided these are on the energy shell 
(p^ = p'^), a minimal requirement for an effective potential. Thus even in the simplest circum¬ 
stances, asking for “the effective electron-electron potential” is too naive. 

More serious is the fact that neither Ucb or Ucg are permissible approximations to an effective 
potential, in the context of a Dirac-spinor description of spin-1/2 particle wave functions. This 
can be most easily seen by including a third particle “3”. The use of a sum such as 17(1,2) -|- 
17(1,3) -|- 17(2,3) to describe the interaction is disastrous from the start, because of the problem 
of continuum dissolution [4,5]. Analysis shows that a theoretically well-founded choice for the 
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second-order potential is neither Ucg nor Ucb, but either of the two operators Vcg or Vgb, 
dehned by 

Vcg = ^++Ugg^++, Vgb = ^++Ugb^++, (14) 

where A++ = A_|_(1)A_|_(2) is the product of positive-energy Casimir-type projection operators for 
the leptons [5,6]. Since the spinors in (8) are eigenfunctions of the A_|_(i), these operators also 
reproduce the lowest-order amplitude. 

All this is by now fairly well known. But a related question has received relatively little 
attention: What is an effective potential which describes the scattering amplitude correctly to 
fourth order? Before considering this, it is useful to study two simpler but instructive cases: two 
spin-0 particles, in scalar QED, and a spin-1/2 and a spin-0 particle, in spinor-scalar QED. 

2.2 Two spin-0 particles 

On using Feynman gauge to write down the one-photon exchange amplitude, one gets 

= 6162(^1 -h Pi) ■ {p2 + P 2 ) A = 6162(20 -h t)/t, (15) 

where a = s—m\—m\. Fourier transformation of the c.m.value of (2.10) yields a term proportional 
to Ugi with an energy-dependent coefficient, plus a contact term proportional to 5(r). Such a 
potential is not suitable for use in a Schrodinger type of equation. In second-order perturbation 
theory it would lead to an ultraviolet (UV) divergence. A potential which is iterable can be 
obtained by first writing in a different form (which does not change its value on the mass 
shell) and then Ending an equivalent operator in r-space, now involving derivative operators [7]. 
One is thereby led to what can be termed a Feynman-gauge-inspired (FGI) potential , which 
in the c.m. system (pA —» Pop,P 2 ^ —^ ~Pop) is 

V[^^^ = z'^Ugz'^ + |/op(PopUcPop/2mAmB)yop, (16) 

with = (1 + Pop/2Ei'’E2^)^'^^ and Pp = { 77117712 /. The corresponding Coulomb-gauge 
inspired (CGI) one-photon exchange potential is given by [8], 

= Vg + Vt (17) 

where, with a curly bracket denoting an anticommutator, 

Vg = yop[{Ey, {FA, UG}}]yop/4:77ii77i2 (18) 


is a relativistic version of the Coulomb interaction and 

Vt = -(l/2)|/op{pA, (pA, (5ij + fifj)Uc}}]yop/4mim2 (19) 

is the potential arising from the exchange of a transverse photon. 
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2.3 One-photon exchange potential in spinor-scalar QED 

To complete the one-photon story, consider the exchange of a photon when “1” has spin-0 and 
“2” has spin-1/2. The Feynman amplitude is then given by 

= -eie2'U/72 ■ (pi + 'p'i)ult (20) 

The corresponding FGI potential is [9] 

vr°‘ = 9rAr(2)C/WAr(2)!,r In?’ ^ (21) 

where 

[/(2) = - a ■ PopUc - Uca ■ Pop)/2mi. (22) 

This is an analogue of for two spin-1/2 particles. There is a corresponding CGI potential 

which I do not write down. Use will be made of (21) in Sec. III. 

2.4 Remarks on the orbit-orbit interaction 

The difference between the two choices, and is connected with the form of the so-called 

orbit-orbit interaction. To see this, note that in the n.r. limit (16) yields as the leading correction 
to Uc an orbit-orbit interaction Uq-o of the form 

= {p^, (5ijUc}}/4mAmB, (23) 

whereas (17) yields 

= (l/2){pr, {p/P, (5y + bfj)Uc}}/4mAmB. (24) 

Note that (24) is a manifestly hermitian form of the orbit-orbit interaction familiar from atomic 
physics, usually described in texts as resulting from the reduction of the Breit operator (13) to 
n.r. form. This is unfortunate from a pedagogical point of view, since the Breit operator refers 
only to spin-1/2 particles. Spin has nothing to do with it! I will return to the difference between 
(23) and (24) shortly. 


3 Scattering potentials: Two-photon exchange 

A key point in the dehnition of any potential, often overlooked, is that one must specify in advance 
how it is to be used. Let us work in the c.m. system of the particles, scattering with total energy E. 
With ho = ho(l) -l- ho(2) the free Hamiltonian for the particles, we will require that the transition 
amplitude computed from 


Top{E) - Veff Veff{E - ho - Veff + le) Ve// (25) 

reproduce the scattering amplitude obtained from field theory to a given accuracy. I consider first 
the case of two spin-0 particles, then the mixed case of a spin-0 and a spin-1 particle, and hnally 
the case of two spin-1/2 particles. 
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3.1 Two spin-0 particles 

In scalar QED there are five Feynman diagrams corresponding to two-photon exchange, which 
contribute to the fourth-order amplitude the two-rung ladder (or ’’box”) graph, and the 

two-rung crossed-ladder graph, a pair of ’’single-seagull” graphs, and a ’’double-seagull” graph. 
The sea-gull graphs are ultraviolet (UV) divergent, but these divergences are taken care of by 
renormalization and do not affect the long-distance behavior. More serious is the fact that the 
’’box” and ’’crossed box” graphs are infrared (IR) divergent. Thus it seems at first sight that one 
cannot even begin to talk about a two-photon exchange potential! 

The resolution of the IR problem lies in the recognition that the iteration amplitude Mj arising 
from Vi.y is now also infrared divergent. This divergence is just the field-theoretic counterpart of 
the fact that in NRQM the second (and higher) Born approximations to the scattering of a charged 
particle in a Coulomb field is divergent (for any value of the momentum transfer). It turns out 
that the difference M 2 ^ — Mj is IR convergent, so that the associated potential ¥ 2 ^ is well-defined 
after all. 

I will only state the form of the result at large r (that is, r large compared to the Compton 
wave length of either particle) and at low momenta. On using to compute Mj one finds that 

[ 7 ], 

+ ... (26) 

where, with k = 6162 / 471 , 

02 ^^^ = /2{mim 2 ), = —7k‘^/6nmim2. (27) 


In contrast, use of yields [ 8 ] 




(28) 


where 

C2^^ = 0, = —Ik"^/Qnmim2. (29) 

The difference between the asymptotic forms (25) and (28) can be traced back the difference in 
the associated forms of the orbit-orbit interactions mentioned above. Thus we see that in the case 
of two charged particles the leading asymptotic behavior of depends on the precise definition 
of Viry. This observation resolves a long-standing puzzle concerning conflicting results for the value 
of 62 . Further, as was noted some time ago by L. Spruch, is classical in character, i.e. if h 
and c are restored, 62 turns out to be independent of h . One should therefore try to understand 
the source of this term from classical electrodynamics. It turns out that this is indeed possible by 
a reexamination of the work of Darwin [1], but I will not enter into the details here [10]. 


3.2 A spin-1/2 and a spin-0 particle 

A similar analysis can be carried out for the mixed case. If one uses to compute the iteration 
amplitude one finds that the spin- independent part of V 27 is 

Rs/*' = k‘^[2^{mi + m 2 )]■V■^ (30) 
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at large r and low energies, which coincides with that for two spin-zero particles. The spin- 
dependent part is essentially a spin-other-orbit type of interaction [9] 

^7°' = -(eie2/47r)^(a ■ £/4m2)[(3mi hm2)/mi{mi + m-2)]r~^ (31) 

The two-photon-exchange spin-orbit interaction therefore decreases as or one power more 
rapidly than that arising from one-photon exchange. Applications of this result to a variety of 
exotic systems, such as pionium and muonic helium, are considered in Ref. [9]. 

3.3 Two spin-1/2 particles 

The computation of for two spin-1/2 particles, even of its asymptotic form for large r, turns 
out to be remarkably complex and has not yet been completed. However, some aspects are known 
and others can be guessed. For example, one can show that if is used for the lowest-order 

potential, then the spin-independent part of V 27 coincides, for large r , with that for two spin -0 
particles and the spin-other-orbit part coincides with that for a spin -0 particle “1” and a spin- 1/2 
particle “ 2 ”, plus a similar term with “1” and “ 2 ” interchanged. At large r and low energies one 
expects that the spin-spin part ¥2^^ has the form 

= k‘^[A{r)ai ■ 02 + B(r)ai ■ fai ■ <72 ■ f], (32) 

with A{r) = ar~^ and B{r) = br~^. The coefficients a and b , and the exponents p and q , 
remain to be determined. For equal-mass particles, such as two electrons, a = where a' 

is dimensionless and likely to be of order unity. To guess the value of p , note that in the large 
r , low-energy limit the spin-spin part coming from Vi^ has one more inverse power of r than the 
spin- other-orbit part. If the same relation holds for two-photon exchange, then p = 5. 


4 Scattering potentials and bound states 

Let us now consider the possible relevance of all this to the computation of relativistic atomic 
structure, one of Ingvar Lindgren’s keen interests. 

The application of QED to compute such structure to an accuracy beyond order a'^Ry, for sys¬ 
tems more complicated than the classic two-body systems (hydrogen, positronium, and muonium) 
began in the mid-hfties, with calculations of the hne structure of helium to an accuracy of order 
a^Ry [11,12]. These calculations, and subsequent ones of spin-dependent effects of order a'^Ry 
[13], were based on a generalization of the two-body Bethe-Salpeter (BS) equation to include an 
external held. The generalized BS equation leads in a somewhat convoluted way to an equation 
of the form 

(33) 

where 

//++ = <(1) + <(2) + L++Hc(l, 2)T++, (34) 

with = tti ■ Pi -f- j3\m. — Za/q the usual external-held Dirac Hamiltonian and L_|_+ the prod¬ 

uct of the associated external-held positive-energy projection operators. This equation, called the 



(external-field) no-pair Coulomb-ladder equation in [12], can be used as the starting point of a sys¬ 
tematic perturbative approach to the energy levels of He and He-like systems. The corresponding 
equation with Uc replaced by Uq + Ub contains not only all a^Ry corrections to hne-structure, 
but a number of corrections of order a^Ry as well. (For large Z these can be expressed in ana¬ 
lytic form. Recently, Lindgren and his collaborators [14] have developed methods for solving such 
equations numerically with sufficient accuracy to enable them to verify the analytic results.) 

However, the further calculations needed to obtain the remaining a^Ry and still higher-order 
corrections, such as those coming from two-photon exchange, are complicated, messy, and have all 
the appeal of a thorny, black box. They contrast strongly with the beauty of scattering-amplitude 
calculations, where among other things one can advantageously use a covariant propagator for 
the photon. It is an attractive idea to try to understand what part of the level structure can be 
thought of as arising purely from those forces which act when the electrons undergo free scattering. 
To put it another way, note that in NRQM one can describe both the scattering of two electrons 
and their interaction within a bound state with one and the same interaction potential Uc- The 
question then is to what extent one can, within the framework of QED, understand the bound 
states in terms of the same two-body effective potentials which describe the scattering accurately. 
It is clear that in QED there are effects which cannot be so described, e.g. the effects of the 
external potential provided by the nucleus during the exchange of photons between the electrons. 

The availability of a lepton-lepton potential which reproduces the scattering amplitude exactly 
to order ( 6162 )^ would be a signihcant step in the direction of such a goal [15]. An early test of 
such a program would be provided, for example, by solution the “free no-pair ladder equation” [5] 

/i++'F = ET (35) 

where 

h++ = hf(l) + hf( 2 ) + I 4 +. ( 36 ) 

Here = cti ■ Pi -t- /djm -|- A+(i)(—ZQ;/ri)A+(i), and V++ is of the form V++ = A++f/(1, 2)A++, 

with A++ the product of free positive-energy projection operators, as in (14). The operator 
U{1, 2) = must be chosen so that if the external potential is turned off, the one-photon 

and two-photon exchange scattering amplitudes are reproduced. (I note in passing that if one 
uses = UcG rather that Ucb one must already include to get the correct a^Ry level 
shift!) On comparing the eigenvalues of (35) with those obtained from, e.g., the Bethe-Salpeter 
equation approach, one could disentangle those dynamical effects involved in the scattering from 
those which are peculiar to the bound-state situation. For Z not too large one would have an 
initial description in the language of conhguration space which is not only quite accurate, but 
which also could actually be explained to students, without requiring that they become experts 
in the arcana of four-dimensional wave equations and all that. So it may well appeal to Ingvar 
Lindgren: As the photo-op at the end of this symposium showed, the results of his pedagogical 
skill and devotion could only be captured by a wide-angle lens! 
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